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The variation of spray characteristics along a radial distance downstream 
of a pressure-swirl atomizer was measured by laser light-scattering technology. 
An analytical model was developed to predict the variation of spray character- 
istics along the radial distance. A comparison of the predicted and experi- 
mental data showed excellent agreement. Therefore, the spray model proposed, 
although relatively simple. Is correct and can be used, with some expansion 
and modification of the prepared model, to predict more complicated spray 
systems . 


INTRODUCTION 

The characteristics of the spray formed by an atomizer In a liquid-fueled 
combustion device are extremely Important to the performance, stability, and 
pollutant formation of the combustor. The spray characteristics should be 
measured accurately, and a model should be developed to predict the variation 
of spray characteristics, droplet trajectory, spray dispersion, and evapora- 
tion history. It Is also Important to understand the Influence of various 
factors on the measurement of spray characteristics. A comparison between 
predicted and measured spray characteristics provides useful Insight Into 
spray combustion. Such research has a significant effect on engineering 
applications and on the fundamental understanding of spray combustion. 

The authors analyzed the variation of spray characteristics along the 
axial distance downstream of a pressure-swirl atomizer. The effect of spray 
evaporation on the variation of spray characteristics Is discussed In 
reference 1 and was presented by J.S. Chin and O.Y. Zhu at this symposium. 
(Paper, entitled The Interdependence of Spray Characteristics and Evaporation 
History of Fuel Sprays In High Temperature Airflows, was unavailable for 
printing at the time of publication.) The Influence of downstream distance on 
the spray characteristics of pressure-swirl atomizers Is caused by the effects 
of spray dispersion and drop acceleration or deceleration. This condition Is 
analyzed In reference 2. These analyses have been partially substantiated by 
some experimental data. Until now there have been few experimental data of 
the variation of spray characteristics along the radial distance downstream of 
the pressure-swirl atomizer. Lee Oodge published his measurement results as 
shown In figure 1 (ref. 3). He shows that both the Sauter mean diameter SMD 
and drop-size distribution parameter N Increase with radial distance. There 
Is no analysis on this aspect available at the present time. 

We are Interested In comparing the experimental and predicted values of 
the variation of spray characteristics along the radial distance. Therefore, 
we carefully measured this variation In a well-defined experimental condition, 
developed an analytical model to predict the variation, and then compared the 
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experimental and predicted data. These comparisons promote our understanding 
of spray dispersion and the factors Influencing drop-size measurement. After 
we have predicted the simple spray system, we will be able to expand and 
modify the proposed spray model to extend Its capability to more complicated 
spray systems. 

Cao and Chin proposed a flat-fan spray model for the fuel distribution 
downstream of a plain orifice Injector under cross airflow (ref. 4). The 
model was validated by fuel distribution measured by gas analysis and drop 
size. It Is obvious that, for the validation of a given spray model, the 
spray characteristics measurement Is reliably accurate and much simpler than 
the fuel distribution measurement. It Is this author's Intention to validate 
a spray model by measuring the spray characteristics along the radial distance. 
The same spray model has also been validated by measuring the spray character- 
istics along different e angles (as shown In fig. 5). 


TEST APPARATUS, ATOMIZER, AND PARTICLE SIZER 

The test section was a rectangular chamber with transparent windows for 
optical drop-size measurement. The test apparatus Included an air system and 
a fuel system, as shown In figure 2. 

The atomizer tested (fig. 3) was a pressure-swirl atomizer taken from an 
existing aircraft gas-turbine engine. The atomizer had a pilot-fuel flow and 
a main-fuel flow. In this test we used the pilot-fuel passage. Because the 
flow rate Is relatively lower with pilot-fuel flow, the obscuration was In the 
right range. 


A Malvern drop-size analyzer model 2200 was used. The principle of the 
particle sizer Is shown In figure 4. The Instrument was set on the Rosln- 
Rammler distribution mode. The direct results from the Instrument are the 
characteristic diameter B and the drop-size distribution parameter N In 
the R-R distribution which Is expressed by 


Q = 1 



(1) 


where Q Is the fraction of the total volume contained In drops that have a 
diameter less than D. 


The Sauter mean diameter SMD can be obtained from 


SMD 



( 2 ) 


where r Is the gamma function. 


EXPERIMENTAL RESULTS 

During drop-size measurement tests, the laser light beam was first 
adjusted so that It was In the same horizontal plane as the central line of 
the atomizer and that It was perpendicular to this central line. At this 
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position (defined as the y = 0 position), the spray characteristics were 
measured. Then the Malvern particle sizer was put on some size blocks of 
known dimension. This changed the height of the laser beam so that the radial 
distance where the particle sizer measured the drop size was changed. The 
summation of the size block was defined as the radial distance y. The dimen- 
sions of the blocks chosen were 2.5, 3.0, and 10 mm. Thus, the following 
radial distance y could be obtained by different combinations of the size 
blocks: y = 0, 5.5, 10.0, 15.5, 20.0, 25.5, and 30.0 mm. The experimental 

results are shown In figures 6 to 13. From these figures It Is clear that the 
Sauter mean diameter and drop-size distribution parameter Increase with radial 
distance until the edge of the spray has been reached. The experimental 
results obtained by the present authors are quantitatively In good agreement 
with Lee Dodge's results. The obvious explanation of the measured results 
shown In figures 6 to 13 Is that with larger radial distances the possibility 
of the laser light beam meeting the large drops Is greater than It Is with 
smaller radial distances. This shows that when we measure the drop size It Is 
necessary to take a sample that Is representative of the spray formed; how- 
ever, It Is Impossible to define the sample which would fully represent the 
spray. The most reasonable sample can be obtained by measuring the spray at 
the y a 0 position. 

Comparing figure 6 with 8 and 7 with 9 shows that, at a higher airflow 
velocity, the change of spray characteristics along the radial distance Is 
weaker. A comparison of figures 8, 10, and 12 and of 9, 11, and 13 shows 
that, with a lower pressure drop across the nozzle, the change of the spray 
characteristics along the radial distance Is weaker. These results tell us 
that when the spray Is expanded less radially the change of spray characteris- 
tics Is also flatter. The reverse Is also true. 

The experimental data obtained In this way were used to check the pro- 
posed spray model . 


PHYSICAL MODEL 

The spray formed by a pressure-swirl atomizer Is shown In figure 5. The 
physical model was proposed with the following assumptions: 

(1) The atomization process Is completed as soon as the fuel leaves the 
atomizer nozzle; that Is, droplets of different sizes start their movement at 
the nozzle exit. Droplets of different sizes have the same Initial velocity. 
Because It Is suitable to use a cylindrical coordinate system to describe 
droplet motion, the tangential velocity and the radial velocity component can 
be combined and treated as one velocity component. 

(2) The airflow Is uniform. 

(3) The air Is at ambient temperature, so droplet evaporation can be 
neglected. 

(4) The Influence of different droplet velocities on the drop-size meas- 
urements made by the Malvern particle sizer Is neglected. 

The Initial drop-size distribution can be expressed by D 0 (SMD 0 ) and N 0 
with the Initial volume fraction (dQ/dD) 0 as 
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At a downstream distance x, the droplets of diameter D move from the origin 
and then distribute themselves on a circumference of radius R(D) while the 
laser light beam of diameter dj only meets some of these droplets: 


2(<P 0 ~ ‘Pi ) 


/dfl\ ‘^2 ~ *1 

\d D ) 2ir 


where 
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The definitions of q> 2 and <p-| are shown In figure 5. 

y 0 - 5 d j i R ' D > i y 0 * 5 d j 


then 


When 


(4) 


(5) 

( 6 ) 


<p 2 


2 


where y Q Is the radial distance and dj the diameter of the light beam. 

The summation of the liquid fraction of the droplets of various diameters 
passing through the laser light beam Is given by 



(7) 


The minimum diameter D m ^ n that the laser light beam can see Is determined 
by 


R ( D m1n> 



That Is, at the downstream distance x, If R(D) < y 0 - 1/2 dj, this droplet 
will not pass through the laser light beam. The relationship R(D) can be 
determined by the droplet motion equation. The maximum diameter that may 
exist In the spray, can be approximated by 


D max = 3.5(SMD) 


If we take the droplets within the laser light beam as a new spray (a 
sample), the liquid volume fraction of droplets of diameter D In this new 
spray (sample) Is 



Since this Is actually the drop-size distribution of the sample (new spray) 
the laser light beam Is to measure, we may obtain Q‘ as a function of D. 
Then, If we assume that the Rosln-Rammler (R-R) distribution can be used for 
data fitting the sample (new spray), we obtain new SMD* and N' for the 
sample. These values, SMD 1 and N‘, change with the distance y 0 . If we 
delete the prime, we obtain the functional relationships SMD = f(y) and 
N = f(y) . 


DROPLET TRAJECTORY EQUATION 


If we use a coordinate system which Is moving with the airflow, then the 
droplet motion equation Is 


2 

dw _ r p a 
dt D 2 


where 

3 

m droplet mass, (ir/6)p 4 D 
liquid density 

w droplet relative velocity 

Cp droplet drag coefficient, Cp = 15/^Re for Re = p a DW/v fl 
p. density of air 

a 

v- viscosity of air 

a 

Equation (9) can be rewritten as 

d ( W o) _ 45 (p a v a w o ) 1/2 n -3/2/w \ 3/2 

dt = "4 p a V“J 

where W Q Is the droplet Initial relative velocity. 


(9) 


( 10 ) 
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Integrating equation (10) yields 


S-- < Bt 


+ 1) 


-2 


(ID 


where 


/ ui \ 1/2 

n 45 ( p a v a W o) .-3/2 

B = 8 Z ° 


( 12 ) 


We know that, at the downstream distance x, droplet velocity, air veloc- 
ity, and the relative velocity of these two have the following relations: 


and 


V = V + (V - v )n- 
x a v xo a'W 

o 


V = v — 
R Ro W 

o 


(13) 


where V x and Vr are the axial and radial velocity components of the drop- 
let, and V xo and Vr 0 are the Initial values of these components. Thus, we 
obtain the droplet trajectory equations 


and 



From equation (11) we have 


(14) 


f 


W_ . t 

W dt " Bt + 1 
0 


(15) 


From equations (14) and (15) we have the radial position (droplet size 0) at 
downstream distance x: 


R(D) = V 


1 


Ro 


B + t 


(16) 
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Time t can be determined by 


XB - V 


t = 


xo 


* [(XB - V xo ) 2 ♦ 4V a Bx] 1/2 


XO' 

2V B 

a 


(17) 


where B Is determined by equation (12). 

From equations (14) and (15), by substituting x and R * = R(D * ) = 

y Q - (1/2) dj, we obtain 


min 


45 (V, M ,) 
8 p. 


1/2 
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(18) 


CALCULATION METHOD AND RESULTS 

In order to calculate the variation of spray characteristics along the 
radial distance the following parameters must be known: 

(1) Axial distance x 

(2) Pressure drop across the nozzle APj, 

(3) Fuel properties and physical properties of air 

(4) Atomizer spray cone angle 

(5) Initial spray characteristics (SMD) 0 and N 0 

The calculation method used In this paper may be summarized as follows: 

(1) Use the measured SMD and N values at y = 0 and the downstream 
distance x as the approximation of (SMD) 0 and N 0 . We know that there Is 
some difference between (SMD) 0 and SMDy_o, N 0 , and Ny_g (ref. 2). So far 
these are the most reasonable approximate values that can be used. 

(2) Calculate the Initial droplet velocity Vj, = y^Apg/pj.)^ where y 
Is the flow coefficient for the pressure drop. For the atomizer used In the 
present research, the spray angle Is close to 90°; therefore, 


( 19 ) 


(3) Determine the D m ^ n value for a different y 0 value by 
equation (18). We assume D^ax * 3 . 5( SMD) . 

(4) Use the Simpson method to solve equation (8), with the R(D) value 
calculated from equations (16) and (17). 
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(5) Determine the Integration 



From this, we are able to obtain a set of the data of Q' and D. 
(6) Assume that Q' can be fitted by R-R distribution: 



By using the least-square method, we can obtain D and N ' ; then, SMD' and 
N 1 can be calculated. These SMD 1 and N 1 values are plotted against y 
In figures 6 to 13. 

The results of the calculations for different airflow velocities and for 
the nozzle pressure drop are shown In figures 6 to 13. These are compared 
with the corresponding experimental data. It Is obvious that the predictions 
are In excellent agreement with the experimental results. 


CONCLUDING REMARKS 

In the present study, droplet evaporation was neglected (discussed by 
J.S. Chin and J.Y. Zhu In the unpublished paper cited previously), and It Is 
not difficult to take Into consideration the effect of droplet evaporation. 

After this model has been validated, It Is relatively simple to predict 
the fuel distribution downstream of a pressure-swirl atomizer (ref. 6). 

One Interesting point that was noticed from the experiments and the anal 
ysls was that the larger droplets are always at the spray edge. Because the 
Ignition spark plug Is also always positioned close to the spray edge, the 
spark might meet the larger droplets, which would harm the Ignition. 

We measured the variation of spray characteristics along the radial dis- 
tance was by laser light-scattering technology. And we found that the predlc 
tlon based on the proposed spray model was In excellent agreement with the 
experimental data. 
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Figure 2. - Experimental system. Plenum chamber, 1; air 
temperature, 2; air flowmeter, 3; test section, 4; exhaust 
pipe, 5; fuel tank, 6; filter, 7; pump, 8; return valve, 9; 
valve, 10; pilot-fuel flow valve, 11; pressure gauge, 12; 
atomizer, 13; main-fuel flow valve, 14; pressure gauge, 15. 





Figure 3. - Atomizer. Pilot-fuel 
flow, 1; main-fuel flow, 2. 
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Figure 4. - Laser drop particle sizer. Laser, 1; beam expander, 2; test 
section, 3; fuel nozzle, 4; lens, 5; photo detector, 6; solid valve, 7; 
amplifier, 8; A/0 converter, 9; data transfer, 10; Interface, 11; 
computer, 12; printer, 13; micro tune, 14. 













Figure 6. - Variation of Sauter mean diameter (SMD) along radial 
distance (y) for x = 70 mm, V a = 9 m/sec, and AP = 12 kg/cm?. 



Figure 7. - Variation of drop-size distribution parameter (N) 
along radial distance (y) for x = 70 mm, V a = 9 m/sec, and 
AP = 12 kg/cm2. 
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Figure 10. - Variation of Sauter mean diameter (SMD) along radial 
distance (y) for x = 70 mm, V a = 6 m/sec, and AP = 8 kg/cm2. 
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Figure 11. - Variation of drop-size distribution parameter (N) 
along radial distance (y) for x = 70 mm, V a = 6 m/sec, and 
AP = 8 kg/cm2. 
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Figure 12. - Variation of Sauter mean diameter (SMD) along radial 
distance (y) for x = 70 mm, V a = 6 m/sec, and AP = 4 kg/cm 2 . 
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Figure 13. - Variation of drop-size distribution parameter (N) 
along radial distance (y) for x = 70 mm, V a = 6 m/sec, and 
AP = 4 kg/cm 2 . 



